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Abstract 



In this contribution we would like to revisit the problems of ferromag- 
netism (F) and antiferromagnetism (AF) in the pure itinerant model. These 
methods can be extended later to the superconducting materials. In our 
model we assume the extended Hubbard Hamiltonian. Transition from the 
paramagnetic state to the ordered state of magnetic nature is decided by the 
competition between kinetic and potential energy in which there is an increase 
in the kinetic energy moderated by the inter-site interactions, and a decrease 
in the potential energy. The competition between these two energies results in 
the existence of critical values of interactions for creating magnetic alignment. 
Only when existing in a given material interaction exceeds the critical value 
for a given type of ordering we can have the alignment of this type. 

The influence of inter-site correlation on F and AF in the presence of 
Coulomb on-site correlation is investigated. The well known Landau free en- 
ergy expansion is corrected for both F and AF ordering to reflect the presence 
of critical interaction in the second order term, which leads to the magnetic 
alignment. 

The numerical results show that the inter-site interactions favor F at the 
end of the band and AF at the half- filled point. In both cases of F and 
AF ordering, these interactions lower substantially the Curie and Neel's tem- 
perature towards experimental data. This helps to remove the paradox in 
magnetism that has persisted for a long time. 
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1. INTRODUCTION 

In the solids there are many phenomena such as band magnetism, metal-insulator transi- 
tion, high temperature superconductivity, which are closely related to the Coulomb correla- 
tion. To describe this correlation scientists have developed what is known today as the single 
band Hubbard model [1-3]. In the original version of this model only the on-site Coulomb 
repulsion, U — |l/r| was included. The use of the simple mean-field approximation 
for this model led back to the Stoner model [4-6], in which the majority and minority spin 
bands are shifted with respect to each other for the quantity AE = m ■ U, where m is the 
magnetization in Bohr's magnetons, and the on-site Coulomb constant U plays the role of 
the exchange interaction. The Coulomb constant U resulting from the Stoner condition for 
creating ferromagnetism is large, i.e. of the order of the bandwidth. On one hand it can 
be justified by the existing strong Coulomb interaction, but on the other hand for such a 
strong interaction, one can not use the mean-field approximation. 

These limitations of the mean-field approximation started the search for a more realistic 
ferromagnetic solution using higher order approximations, e.g. the Hubbard I approximation 
[1], coherent potential approximation (CPA) [7,8]. While the Hubbard I approximation 
describes correctly the atomic limit (t = 0) it is not correct for the intermediate U/t. For 
these values the coherent potential approximation (CPA), which is equivalent to Hubbard 
III decoupling of the Green functions, is better. Unfortunately, after a thorough analysis it 
was concluded that the Coulomb repulsion U, no matter how strong, does not lead to the 
ferromagnetic order (see e.g., [9,10]). The reason is that it does not yield the spin dependent 
band shift or spin dependent band narrowing necessary for a ferromagnetic ordering. 

Both Hubbard I and CPA approximation have two types of particles; electrons moving 
between empty sites and electrons moving between sites occupied by electrons of opposite 
spin, separated from each other by the energy U and forming lower and upper Hubbard's sub- 
bands. This is the two-pole approximation in the language of Green functions. Base on this 
scheme is the newer model called "spectral density approach" (SDA) [11,12]. The magnetic 
phase diagrams calculated by this method are more realistic and the values of the Curie 
temperatures also make sense. The main advantage of SDA is to obtain a spin-dependent 
band shift necessary to obtain the ferromagnetic ordering. On the other side the SDA method 
is the linear combination of two 5 functions leading to the real self-energy, and lacking the 
quasi-particle damping. Therefore, Nolting and co-workers have proposed a combination of 
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SDA and CPA called "modified alloy analogy" (MAA) [13]. This approximation has brought 
the self-consistent ferromagnetic solutions, and the complex values of self-energy. The defect 
of this method is the small range of concentrations for which one has the ferromagnetism 
and also a slightly unclear derivation of the method. 

Another approach to solving the Hubbard model is computer simulation. The initial 
attempts of the Monte Carlo simulation applied to the basic Hubbard model also did not 
give the ferromagnetic ground state [14]. More recently, the new dynamical mean- field 
theory (DMFT) [15,16] has been developed for a direct computational simulation of systems 
with correlated electrons on a crystal lattice. This method has the exact solution in the non- 
trivial limit of an infinite coordination number [17]. The results have been obtained by using 
quantum Monte-Carlo (QMC) simulation [18-20] and the mean-field Green function theory. 
Use of this dynamical mean field theory (DMFT method), has introduced a significant 
progress in the theory of ferromagnetism. The results (see e.g. [21,22]) show the existence of 
ferromagnetism but at much lower temperatures than those coming from the Hartree-Fock 
approximation. Such results could remove the problem known as a "magnetic paradox" , i.e. 
the Curie temperature coming from the interaction constant (e.g. U) - fitted to obtain the 
correct magnetic moment at T = K - as being much too high. The later attempts along 
the computational lines (see [16,20,22-25]) are still inconclusive, mainly for the reason of 
the limited size of the computational systems which seems to be too small too represent the 
bulk materials. 

In recent years there has been a very fruitful analytical development in describing the 
magnetism by extending the Hubbard model and including into it the band degeneration 
and the inter-site Coulomb interactions. The influence of inter-site interactions on the band 
magnetism has been studied by many authors [26-37]. They have shown how different inter- 
site interactions would affect the ferromagnetism in the presence of the on-site Coulomb 
repulsion. This also will be the subject of the present chapter. 

On the experimental front in the last few years there has been discovery of a triplet 
superconductivity in UGe2 [38], URhGe [39] and ZrZn2 [40] materials. Although the transi- 
tion to superconductivity takes place in very low temperatures (~0.8 K), the mechanism of 
this phenomena is very interesting, since the superconducting pairing is (probably) stimu- 
lated by the Hund's interaction (see e.g. [41]), the same one which stimulates weak itinerant 
ferromagnetism. It was also found that the growing magnetic moment would suppress the 
superconductivity. The explanation of these facts would be quite simple; the Hund's inter- 
action stimulates both F and SC ordering, but if the ordering of a given type already takes 
place then it precludes the other type of ordering. 

Another type of new phenomenon which includes the band magnetism is the well docu- 
mented coexistence between superconducting ordering and itinerant antiferromagnetism in 
high-Tc cuprates, e.g. YBa2Cu306+a;, La2^a;(Sr,Ba) x Cu04. Since in these compounds the 
magnetic moment is stronger it suppresses the superconductivity until the moment itself 
disappears (see [42-44]). 

Because of all these reasons of all these theoretical problems and experimental facts we 
would like in this chapter to revisit the problems of ferromagnetism and antiferromagnetism 
in the pure itinerant model. These methods can be extended later to the superconducting 
materials, where the magnetic ordering competes rather than cooperates with superconduc- 
tivity. 
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2. THE MODEL 



The general model of itinerant magnetic ordering is based on the extended Hubbard 
quasi-single band Hamiltonian of the following form 



(1) 



H = - Y. tij ( c ta c ja + h - c ) - A*o E hi - F Y, n ia n ir7 + U £ h^h il + V £ h i h j 

<ij>cr i ia i <ij> 

+J £ ctcpd a/ c j(T + J' E (cJcJcjiCjt + /i.e.) 

<ij>a,a' <ij> v 



where /xo is the chemical potential, (c^) creates (annihilates) the electron with spin a on 
the i-th lattice site, hi a = cf a Cicr is the particle number operator for electrons with spin a 
on the i-th lattice site, hi = h ia + hi_ a is the charge operator, n ia is the average number of 
electrons on sites i with spin a. 

In this Hamiltonian we included not only the dominant on-site Coulomb correlation, U, 
but also the on-site Hund's field F. Such a field can exist only as the interaction between 
different orbitals in a multi-orbital model. We assume the single band, but composed from 
identical orbitals, which are fully degenerate i.e. have the same density of states and the 
same electron occupation (see [10]). In such a band the effective exchange field can be 
expressed as; F = (d — 1) • J in , where J in is the exchange interaction between different 
orbitals within the same atomic site, and d is the number of orbitals within the band. As 
a result our model is a quasi-single band model. The intra-atomic Hund field in Eq. (1) is 
already expressed in the Hartree-Fock approximation, which will be justified only for small 
values of this interaction. In addition we have three explicit inter-site interactions [1,45,46]; 
J-exchange interaction, J'-pair hopping interaction, V-density-density interaction. To avoid 
a large number of free parameters we will assume later on in the numerical analysis that 
J' — J, V — 0, which will leave us with only two parameters; J and parameter S representing 
kinetic interactions (see Eq. (5)). 

The spin dependent correlation hopping tfj depends on the occupation of sites i andj, 
and in the operator form can be expressed as 

t£ = t(l - rV CT )(l - nj- a ) + h [hi- a (\ - hj_ a ) + nj- ff (l - hi- a )\ + t 2 n i - ff n j - ff , (2) 

where t is the hopping amplitude for an electron of spin a when both sites iand j are empty. 
Parameters; t\, t%, are the hopping amplitudes for an electron of spin a when one or both of 
the sites % or j are occupied by an electron with opposite spin, respectively. Quite recently, 
several authors have suggested that the expected relation; t > t\ > t 2 , may be reversed for 
large enough inter-atomic distances; t < ti < t 2 (see Refs. [47] and [48]). This concept would 
fit the results of Gunnarsson and Christensen [49], who for the heavier elements (e.g. 3d or 
4f) claim growing hopping integrals with increasing occupation. 

Including the occupationally dependent hopping given by Eq. (2) into the Hamiltonian 
(1) we obtain the following result 



H = - £ [t - At (hi_ a + n j _ CT ) + 2t ex h i _ a h j _ a \ (cl a c ja + h.c.) - /j, E hi ~ F £ n^n^ 

<ij>a i ia 

+U £ hi^hii + V £ hihj + J £ ct a c+ a ,c i(jl c j(J + J' £ (c^c^Cjic^ + h.c.J 

i <ij> <ij>a,a' ' <ij> 

(3) 
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where 

t + t 2 

At = t-h , t ex = —^--t 1 . (4) 

In this form it is quite visible that the kinetic interactions: the hopping interaction, At, 
and the exchange hopping interaction, t ex , are also the inter-site interactions. 

We assume that ti/t = S and £2/^1 = S±. In general these parameters are different and 
they both fulfill the condition S < 1 and Si < 1 which is equivalent to t > ti > t 2 (see Ref. 
[47]). For simplicity we will assume now additionally that S = Si. With this assumption 
we have the relationships 

At = t-ti=t(l-S) , t ex = t -^-ti= t -(l-S) 2 . (5) 

In the Hamiltonian (3) there are many inter-site interactions: At, t ex , J, J', V, for which 
we will use the modified Hartree-Fock (H-F) approximation. For the two- and three-body 
terms of Hamiltonian we will adopt the procedure introduced by Foglio and Falicov [50], 
Aligia and co-workers [51], and Hirsch [27]. Neglecting superconductive averages of the type: 
c i- c jo-')> ( c icrCj a '}, and the spin-flip terms (c^Cj-^ we obtain 

- e.g. for two-body term 

hihj = (n iT + h^) (n n + n n ) (n, T + h^) (n n + n n ) + (n iT + n n ) (n n + n n ) 

- E (( C ta C jv) Cja C ia + h - C -) + COnst ' 

- e.g. for three-body term 

c^c n h tl n jl = c+c^cf^c^c^ cJc jT ((n^) (h n ) - (cjc^) (c^c^)) 

-4l c jl ( c iW) { c ti c n) ~ c n c n ( C *W) { c ti c ii) • (7) 

+n n (c+ T c jT ) (n jir ) + n n (cJc jT ) (n^) + const 

The main point of this approximation is to retain the following inter-site averages: I a = 
cfpCja^, in addition to the usual on site averages: nj CT = (fii a ), which contribute to the 
Stoner field. 

After performing this modified Hartree-Fock approximation on all inter-site interactions: 
At, t ex , J, J', V, we obtain the following simplified Hamiltonian 

H = - E %U { C taCja + h.c.) - 110 E + E M >™ + U H "iT*M > (8) 

<ij>cr i ia i 

where t^- = t ■ b a is the effective hopping integral, with 6°given by 

b a = 1 - {n^ a + rij-a) + 2 f -j- (ni- a nj- a - l\ - 2I a I^ - - — I a - - + - I_ a , 

(9) 

and M? is the spin-dependent modified molecular field for electrons with spin a expressed 
as 
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Mf = -Fn ia - J rij* + V E + + 2*At/_ ff - 2t ex I_ a ]T n j<r , ( 10) 

3 3 3 

where z is the number of the nearest neighbors, J2j is the sum over nearest neighbors of the 
lattice site %. The above molecular field is the sum of the on-site contribution F and the 
inter-site contributions, which will be different for F and AF. 

2.1 INTER-SITE AVERAGES I a = (cf a c j(! 

Parameter I a can be calculated from the average kinetic energy. In the case of ferromag- 
netism one can write the following expression for the average kinetic energy of +a electrons 

(l< a ) 



( K<7 ) = -tiff E = -Kffla = -DZffl. , (ii) 

<ij> 

where D"^ = zt u e ^ is the half band-width of the +a electrons. 
The kinetic energy can be calculated straightforward as 

D 

(TO = J f [e a (e )} ■ e° (f) ■ p° (e°) ds° = -D^ , (12) 

with e a (e°) = b a e° being the energy of the deformed band, e° is the energy of unperturbed 
band with all inter-site interactions being equal to zero, and f \e a (e )] = — ; — \ — —, — is the 

1+e ' 

Fermi function. Eq. (12) means that we can calculate the average product of two nearest 
neighbor operators: I a — — (K a ) /-De// • 

In the case of zero temperature, T = K, and constant density of state (DOS): p° (e°) = 
const = jj for — D < e° < D (D = zt - unperturbed half-bandwidth), we obtain from Eq. 
(12) that 

I a = n a (1 - n c ) . (13) 



This approximation suggests that we can treat the average I a = (c^Cj^J as the probabil- 
ity of electron with spin o hopping from the % to the j lattice site and back. More precisely, 
it is given by the average of two products. One is the product of probabilities that there 

is an electron with spin a on the i site and that the j site has empty states: -""^"-^ H -' T \ , 

and the second one is the probability of the opposite jump; ^q^~f z^y - The quantity hi a 
is the average number of electrons in the sub-band o on sites % and n* CT is the total capacity 
of the sub-band a on site %. The total hopping probability, called I a , is given by the average 
of above probabilities 



nJ^jg - n ja ) h ja {n\ a - h ia ) 
n ia + {n) a - h ja ) h ja + (h\ a - n ia ) 



(14) 



For the weak correlation (U « D), when the band is not split by the Coulomb repulsion, 



we have n\ a = 1 and n ia 
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Ia = 



2 



n ir7 (l-n ja ) n ja (l-n ia ) 



n ia + (1 - n ja ) n ja + (1 - n ia ) 



(15) 



Expressions (14) and (15) give a general form of quantity I a , which depend on the type 
of magnetic ordering and on the strength of Coulomb on-site interaction through the mean 
values n ia and fij a . 

At first let us consider the quantity l a in the ferromagnetic state when one has the 
obvious relation; h ia = fij a . In the case of a weak correlation (U « D) the band is not 
split and for the electrons with spin a we have; h ia = n ia = n CT , and expression (15) can be 
simplified to Eq. (13). 

In the case of strong correlation (U » D) when the band is split into the lower and the 
upper Hubbard sub-band, quantities hi a and n\ c , will depend on the location of the chemical 
potential. When n < 1 the chemical potential is located in the lower Hubbard sub-band and 
we have; h ia = n a , n\ a = 1 — n_ a , what gives for I a the result 

1 - n_ a 

-for the upper Hubbard sub-band (n > 1) we have: n ia = n a — (1 — n_ CT ) —n—1, n\ a = n- a 
and I a is equal to 

T _ ( n - !) (! ~ w <t) ( 17 ) 



In the antiferromagnetic state the crystal lattice will be divided into two interpene- 
trating sub-lattices; a, f3 with opposite spins, and with the average electron numbers equal 
to 

n±m g n=Fm 

n± a = n± a = — - — , n'± a = n T<T = — - — , (18) 

where m is the antiferromagnetic moment per atom in Bohr's magnetons. 

The indices i,j belong to the neighboring sub-lattices a,f3 with opposite magnetic mo- 
ments. In the result for AF ground state the parameter; l a = 7_ CT = Iaf is spin independent 
(i.e. it does not depend on the first power of antiferromagnetic moment m) 

The parameter Iaf for the weak correlation is calculated from Eq. (15) using values of 
n a{fj) from Eq> ( lg ) 

2n — n 2 — m 2 , , 

Iaf = — jz 2T" • 19 

4(1 -m 2 ) 

For the strong correlation and the chemical potential located in the lower sub-band 
(n < 1), we obtain from Eq. (14) that 

(2n-n 2 -m 2 )(l-n) 

J-AF = 77^ ^ 5 ' ( 2() ) 

(2 — n) 2 — m 2 

and for the upper sub-band 
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(n - 1) (2n - n 2 - m 2 ) 
-Uf o n 



n 2 — m 2 



(21) 



After calculating the quantity I a for F and Iaf for AF ordering we can find the bandwidth 
parameter, b a , and the modified molecular field, M£ , which are present in the Hamiltonian 
ofEq. (8). 

2.2 BAND-WIDTH, MOLECULAR FIELD AND ELECTRON OCCUPA- 
TION 



2.2.1 FERROMAGNETIC STATE 



In the ferromagnetic state the electron concentration is the same on each lattice site: 
nia = fija = n a , and the spin dependent parameter I a in the case of weak Coulomb cor- 
relation is given by Eq. (13), and Eqs (16) or (17) in the case of the strong correlation. 
Inserting those values into Eqs (9) and (10) we have 

b° = l-2^n^ + 2 t f(n\-I 2 a -2IJ^)-^I a -^I^ , (22) 



M? = M IT = -(F + zJ)n < r + zVn + 2zI- <T (At-t ex n ff ) . (23) 
Using Eq. (5) we can simplify Eq. (22) to the form 

^ = l_2(l-5)n_ CT + (l-5) 2 (n 2 _ CT -/ 2 CT -2/ CT /_ CT )-^^/ CT -^±^/_ CT . (24) 

The effective total field in a ferromagnetic state, F t ^ t ,can be found now from the Weiss 
assumption that the energy shift between both spin sub-bands, AE, is equal to this field 
multiplied by the existing magnetic moment 

AE = M~ a - M a = Ff ot ■ m , (25) 

which after assuming the Eq. (13) for I a gives 

2 2 
777, 

F* a = F + z[J + t ex + 2At(l-7i)] . (26) 

In analyzing the Hamiltonian (8) for the appearance of ferromagnetism we will use the 
level of approximation which will modify the DOS, in addition to its shift by the molecular 
field. 

After Fourier transform of the kinetic energy the Hamiltonian (8) takes on the following 
form 

H = E ( £ l + ^ - ri^ka + fHafli-a (27) 

with the spin dependent electron dispersion relation given by 
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% = eW , (28) 

where b a and M CT for the ferromagnetic state are given by Eqs (22) and (23), and e° k is the 
initial dispersion energy of the electron (without the inter-site interactions) 

e° k = -tJ2 e* k ( R i- R j) . (29) 

<i,j> 

The Hamiltonian (27) is solved in the CPA approximation (see [8,52]), which produces 
the following equation for the on-site self-energy S CT (e) 

(1 - i + tk(?) + 1 - (u- hm = • m 

with the spin dependent Slater-Koster function F a (e) which can be written as 

F '"»4?r^jiF^ ' (31) 

This function F a {e) can be expressed by the unperturbed function 

• (32) 

by the help of the following relation 

F A s) = L F J ^±^Jk \ , (33) 



which becomes the standard CPA relation: F a (e) = F (e — M CT + /x — E CT ), when all the 
inter-site interactions are zero and as a result b a = 1. The above defined Slater-Koster 
function was obtained from the solution of the CPA equation and depends on the spin 
orientation and on the on-site interaction, U, and the inter-site interactions. 
For the unperturbed Slater-Koster function we have the relation [8] 

p(e) = --lmFo(e) . (34) 

71 



Using the identity 



p a (e) = —ImF a (e) (35) 

7T 



one can calculate the perturbed DOS depending on the on-site Coulomb correlation U 
through the self-energy S CT , and on the inter-site interactions through the band- width pa- 
rameter b a . 

Using this DOS one can write the expression for electron occupation in the ferromagnetic 
state 

°° de 
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where p a (e) is the spin dependent density of states given by Eq. (35), and the previous 
CPA equations. 

The electron occupation number and the magnetization are given by 

n = n a + n_ CT , (37) 

m = n a — n_ a . (38) 

2.2.2 ANTIFERROMAGNETIC STATE, DIAGONALIZATION 

In the antiferromagnetic state the magnetic moment on the nearest lattice sites is 
opposite with respect to each other (see Eqs (18)), the quantity/^ = I_ a = I AF , and the 
bandwidth reduction parameter b a from Eq. (9) is spin independent 

w u-° h^F i At , tex ( n 2 -m 2 2 \ 2J+ J'- V T 

b =b =b = 1 —n + 2— I M AF \ I AF . (39) 

After expressing the kinetic terms by a factor S defined in Eq. (5) it takes on the form 

b AF - 1 ~ (1 ~ S)n + (1 - Sf (r^g- - 3/l F ) - 2J+ f " V Ia F • (40) 

The generalized (modified) molecular field from Eq. (10) depends on the spin and on 
the sub-lattice index: a or (5. For the sub-lattice a we have 

M a a = -Fn a ~z(J + 2t ex I AF ) n_ CT + zVn + 2zMI AF , (41) 

and for the sub-lattice (3 

= -Fn_ a -z(J + 2t ex I AF ) n a + zVn + 2zAtI AF 

M$ = M-° ■ [ > 

Now, we can calculate the effective total field on the sites a, (3 in a similar way to the 
case of ferromagnetism 

F , != MZ^k = ^r =F . z(J+2tejAF) . (43) 

m m 

This equation shows that the positive inter-site exchange interaction, J, and exchange- 
hopping interaction, t ex , are opposing AF. It is contrary to the case of ferromagnetism (see 
Eq. (26)), where both these interactions (when positive) are helping the ferromagnetism. 

As mentioned above we assume the type of the crystal lattice, which can be divided into 
two interpenetrating sub-lattices: a, j3, with the average electron numbers given by Eqs (18). 

For the AF ground state we use the diagonalization of Plischke & Mattis [53], Brouers 
[54] and Mizia [55]. The Hamiltonian (8), after taking into account Eqs (39)-(42), will take 
on the following form 
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H = -teff E ( a £P* + Pi* a i°) + E - to + UnlJ) hi , (44) 

(7 = 

where a^.(aj (T ) and P^(Pi a ) are the creation (annihilation) operators for an electron of spin 
a on the sub-lattice a and (5 respectively, h] a = 7^7^ is the electron number operator for 
electrons with spin a on the sub-lattice 7 = a, (3, t e ff — t ■ b AF is the effective hopping 
integral. 

Now we will derive the equations for dispersion relation, particle number and magnetiza- 
tion in the antiferromagnetic state using Hamiltonian (44) and the Green function technique. 
For the on-site interaction, U, we use the standard CPA approximation, which for simplicity, 
in further analysis, is treated in the weak and strong correlation limits. For the weak corre- 
lation we use a first order approximation in interaction constant over the bandwidth. This 
is equivalent to the Hartree-Fock approximation. The second case is the high correlation 
approximation, U » D, , which will be easy to extend for the arbitrary strength of the 
on-site interaction later on. 

The main idea of the CPA formalism [8] is used now. We split the above stochastic 
Hamiltonian (44) into a homogeneous part 

#o = -teff E ( a tPiv + - f E « + hi) + E Khl + E ^ht , (45) 

iu iu iu iu 

and a stochastic part 

^ = E(C-^)C + E(Kf-^)nf CT , (46) 

iu iu 

where 



/i = fx - zVn - 2zAtI AF (47) 

is the effective chemical potential, are the self energies on sites 7 = a, (3 for electrons 
with spin a and V"^ are the stochastic potentials given by 

^a(B) \ K W) = -Fn%W -z(J + 2t ex I AF )n^ . + , , , ... + . P?f ] = 1 - n a _f 

(48) 

The self-energies fulfill the CPA equations 



2 T/7 _ V7 

E^X \ = 7 = «(/?), (49) 

with the Slater-Koster function F^(e) for the sub-lattice 7 in the following form 

*7(e) = ^EG27M) > (so) 
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and G^ 7 (e, k) given by Eq. (55) below. 

After dropping the paramagnetic part the stochastic potential (48) takes on the following 
form 



where 



ya(fS) = 



n ±m 



= Ut a 



a(J3) 



, with 



Pi 



a(J3) 



1-n 



2a 



r 2a — ' i =F0- 



7 AF 
tot 



rn 



, and F t A ot F = F - z (J + 2t ex I AF ) . 



(51) 



(52) 



This simplification of the Eq. (49) allows to see that the self-energies fulfill the relations; 
E" = E^ CT = E + , and S° cr = E^ = E_. Therefore the solution of the problem can be 
reduced to finding only the self-energies, S±. The interpretation of the parameter A as the 
antiferromagnetic energy gap will be given below. 

We transform Hamiltonian (45) into the momentum space, and use it in the equations 
of motion for the Green functions 



e((A;B)) e = {[A,B] + ) + {{[A,H ]_;B)) £ , (53) 
where (A,B) G (a^, a ka , (3 ka , Pkaj- As a result the following equations are obtained 

(54) 



e + fi - E° -e k 
-e k e + n - E^ 



G?(e,k) Gf(e,k) 
G**(e,k) GW(e,k) 





1 o" 




1 



with e.g. Gf{e, k) =« a ka ; a+ », G**(e, k) =« (3 ka - a+ ». 

Solving this set of equations we arrive at the following expressions for the Green functions 



e + a- E^ (a) 

[G (e eff , k)-G (-e eff , k)] 



with 



2\ £ + n-^ 



e eff = v /(e + /i-E + )(e + /i-E_ 



and G(e e ff, k) given by 



G (e eff , k) 



£eff - £k 



Remembering that 



pi (e) = — Im [Fl (e)} 

71 



we obtain 



Pa 



m = _il m , 



£ + jl - S 



(3(a) 



( ^eflf \ 



tt W^-E? (/3) b AF " UAFj- n 



ImF 



± 



(55) 



(56) 



(57) 



(58) 



(59) 
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where F (e) is the unperturbed Slater-Koster function given by Eq. (32). 
Using Eq. (50) we obtain the expressions for electron numbers 

oo oo 

"J=/TT^(4) Im K( £ »' d£ = ^/TT^(4) Im i G "( £ ' k »l d£ ■ < 60 > 

— oo k — oo 

Treating the self-energies as the solutions of Eq. (49), and using it only in the first 
order approximation with respect to the Coulomb on-site repulsion, U, we obtain 

2 ~ n 
Z« a = X ± ^J2P« ±a Vf ±a = [U-F-z(J + 2t ex I AF )]-TA , (61) 
i=i z 

where the antiferromagnetic energy gap A is given by the expression 



A = 



m 



U + F t A J\^ . (62) 



This expression is different from Eq. (52) by adding to the effective field the Coulomb 
repulsion, U, which in the case of the weak interaction is also treated in the Hartree-Fock 
approximation. 

Neglecting the paramagnetic part in self-energies we have 

E ± = =fA. (63) 
Inserting these self-energies into Eq. (55) and redefining e + p as e we obtain 



GT m (e,k) = y £ -^[G(e eff ,k)-G(-e eff ,k)} . (64) 

From this equation using Eqs (50) and (58) one can obtain the following relation for the 
DOS on sites a, f3 



P^He) = f^p(e eff )= P± (e) . (65) 

In this form it is evident that the AF DOS will vanish between energies —A and +A 
(zero energy was assumed in the center of the atomic band or approximately at the atomic 
level). Hence, the quantity A is the AF energy gap. The schematic shape of the DOS is 
shown in Fig. 1. 
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FIGURES 




One can be prove that if the old distribution is normalized to unity; 



D 



J p(e°)ds° = 1 , (e 



-D 



then the new one given by Eq. (65) is also normalized. 
Since p+(e) = p_(— e), one can write that 

J p + {e)de= J -\p+(e)+p-(-e)]de = - J [p+{e) + p_(e)] de 
where D e ff = zt ■ b AF . Inserting 



P+{e) + P-(e) = I r + \ p^ir ) P (^//) = ^—p (e e// ) (e 



e + A V e- A 
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and using the identity obtained from Eq. (56); de = ^^de e jf, one obtains from Eq. (67) 

(69) 



D 



j p+(e)de= J p(e eff )de 



eff 



1 . 



Restricting as before only to the first order approximation for the self-energies (Eq. (63)) 
in Eq. (60) we obtain the following expressions for electron numbers 



^T,[ p k a f(E k )+prn-E k ) 

where f(E k ) is the Fermi function 

f(E k ) = 



(70) 



1 



1 _|_ e (E k -ti)/kT ' 



(71) 



E k = yftT* , ^ = ^4 , (72) 

and 

Pt = ^ (73) 

is the occupation probability of state (k, ±cr). 

The expression for electron occupation in the antiferromagnetic state in the presence of 
external field H will be different from Eq. (36), which was suitable for ferromagnetism. 
Now, the external field will act differently on both spins ±cr in two sub-lattices a and/5. 
Therefore its energy has to be included into the energy of the internal exchange field in the 
process of diagonalization. It will result in the following formula for the electron occupation 
on the sub-lattice a 

oo 

/(]£ 
P° ( £ > i + e (e-,)/kT > where x ° = aA + °VbH- (74) 

— oo 

The precise form of the density of states p"^ (e) comes from Eq. (59) (or (65) in the 
case of the Hartree-Fock approximation). As it can be seen, it depends on the initial DOS 
and on the result of the diagonalization of two sub-lattices Hamiltonian. 



3. FREE ENERGY, STATIC MAGNETIC SUSCEP- 
TIBILITY 

For a simple magnetic material, the free energy can be written 

F m — F + a 2 m 2 + a 4 m 4 + ... = F + /i|m 2 /x + a^ 4 + ... (75) 
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with x being the magnetic susceptibility. 

The electron occupation of the ferromagnetic state in the presence of an external 
magnetic field H is given by the expression (36) with added shift from the magnetic field 

" £ (76) 



1 _|_ e {e-ti+M°-an B H)/kT ' 



where \ib is the Bohr magneton. 

The magnetic susceptibility is given by 

(d{n a -n- c )\ 

*="*H^L. ■ (77) 

As a result of inserting Eq. (76) into it we arrive at the following equation for the static 
ferromagnetic susceptibility 

2f%h .... 
1 — K — F&I T ' (78) 

where K is the correlation factor defined mathematically by the following equation 



TV- r dpa-(e) <fe r dp- a {e) <fe_ 

am, 1 „(e-M+M^)/fcT J a m fe-u+M 

1+6 > (79) 



dm 1+ (e-n+M<r)/kT J dm (e-y,+M-<?) I kT 

W(e) de 

^ J 9m uJe-p+JCI/it 
— oo ±Te 

F^is the total Stoner exchange field, which on the base of Eq. (25) can be written as 



r tot 



dM a 



dm 



0M a 



„ dm 

m^O 



dM° 

= -2 



dm 



i 

m-»0 



(80) 



and 

It= J Pm =o(e)P T (e)de, P T (e) = = f ( e ) e M/fcT -^ , (81) 

— oo 

/ (e) is a Fermi function. 

Returning with the expression (78) to Eq. (75) we arrive at a new form of Landau 
expansion which is 

F m = F + m 2 1 ~ K ~ F * ot It + a 4 m 4 + . . . . (82) 

The critical (minimum) value of the total Stoner field creating ferromagnetism; 
(F^Jmin = F t ^ t , is obtained from the zero of susceptibility denominator (Eq. (78)) or 
zero of the numerator in Eq. (82), around which changes the sign of the second order term. 
Hence, we arrive at the well-known Stoner criterion (see [56,57]) 

FZ = ^ , (83) 
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which is modified here by including the correlation factor K. 

Using the Stoner condition/ 7 ^ = (1 — K)/It, back in the expression for susceptibility, 
Eq. (78), and in Landau expansion (82) we arrive at their corrected forms 



2/i 



X = pcr t? , (84) 

r tot ~ r tot 

2 

F m = F + — (FZ-F tot ) + a 4 m 4 + ... , (85) 

since > 0, the existence of the minimum in the above equation for F m will depend on 
the sign of 02 = \{Ftot ~ Ftot) or equivalently on the sign of the difference: F^ t — F tot . As 
a result, the simple Stoner criterion (Eq. (83)) will lead to the ferromagnetic instability 
with nonzero magnetic moment being given by the minimum of F m curve in Fig. 2. This 
minimum appears only when existing in a given material, Ff ot > F^ t , the F^ t being given 
by condition of Eq. (83). It is not mistake that in Eqs (84) and (85) we have F to t without 
the upper-script F, since as we will see below those two equations are valid without any 
changes for both F and AF ordering. 

To calculate the static antiferromagnetic susceptibility we will apply Eq. (77) to the 
susceptibility of the sub-lattice 1 

X = »b\ { °- a) \ , (86) 



dH 



H=0 



and we will use expression (74) for the electron occupation on sub-lattice a in the pres- 
ence of an external field, arriving at the following equation for the static antiferromagnetic 
susceptibility 



where the correlation factor K x in the antiferromagnetic case can be written as 



X 1 _ u jpAF ' (87) 





dp«(e) dp%(e) \ de J dp« a (e) 

Q x a Q x -a ) 1 + e (e-fJ.)/kT J fao ^ 



(88) 

Using the condition for the zero of susceptibility denominator we arrive at the critical 
value of the total field 

F -*^y • (89) 



x This susceptibility is a susceptibility of one sub-lattice. To find its relation with the total ex- 
perimental susceptibility one has to take into consideration the direction of the external field with 
respect to the easy magnetic axis (see e.g. [58], [59]) 
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After inserting this value into the sub- lattice susceptibility (Eq. (87)) and next to the 
formula for Landau expansion (Eq. (75)) we obtain again the same form of the susceptibility 
as for ferromagnetism (Eq. (84)), and of the free energy expansion as given by Eq. (85) for 
ferromagnetism. Thus the susceptibility of Eq. (84) and the energy expansion of Eq. (85) 
are the same for F and AF. The difference between these two orderings lies in the different 
values of critical (minimum) field. For F this field is given by Eq. (83) and for AF by 
Eq. (89). In both cases when existing in a given materialF to< > F^ t , then the second 
order term in energy expansion vs. magnetization becomes negative and we have a nonzero 
equilibrium value of the ordering parameter m, see Fig. 2. 




FIG. 2. Dependence of Landau free energy on the magnetization m for various values of 



a 2 



per 
r tot 



_F(AF) ™ r 
r tot r tof 



r tot J ■ 



Curve (a) 
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r tot ^ r tot' 



curve (b)- F t F J AF) 
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4 FERROMAGNETISM 

4.1 ONSET OF FERROMAGNETISM 

To compute the ferromagnetic criterion from Eq. (83) we have to calculate the correlation 
factor, K, defined in Eq. (79). We will calculate it now in the CPA approximation, assuming 
that we know the self-energy E CT . The total change of the DOS with magnetization, K, which 
enhances the possibility of creating magnetic ordering, is given from Eq. (79) as 

K = Kjj + K b , (90) 

where 
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v o + F d Pv( £ ) d ^ ft \j /m \ 

^ = 2 7 ^s7^ /(£)& (91) 

— oo 

and 

— oo 

Factor Kjj describes the role of on-site interaction U, and the role of inter-site inter- 
actions in creating magnetization. 

The schematic depiction of the DOS deformed by the on-site correlation Uis shown in 
Fig. 3. 
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FIG. 3. Schematic DOS showing the influence of the strong on-site Coulomb correlation, U. 
The paramagnetic DOS for both spins, ±<r, are solid lines. At U which is strong enough to split 
the band into two sub-bands, lower sub-bands have the capacity of 1 — n _cr for +a electrons, and 
1 — n a for —a electrons. The changes in the spin electron densities integrated over energy are the 
shaded areas in this figure, and they are equal to the correlation factor Kjj. The shift between +cr 
and —a electrons is created by the assumed exchange field. 



The schematic depiction of the DOS deformed by the inter-site correlation is shown in 
Fig. 4. 
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FIG. 4. Schematic DOS showing the influence of the inter-site interactions. The paramagnetic 

DOS for both spins, ±cr, are solid lines. The inter-site interactions change the relative width of 
the bands with respect to each other (described by the b a factors, see Eq. (22)). The Stoner field, 
which would displace the bands with respect to each other, is assumed to be nonzero. The shaded 
areas in this figure are the correlation factor Kf,. 



Comparing Figs 3 and 4 one can see that the on-site correlation causes an increase in 
capacity of the majority spin band (or lower sub-band in the case of split band), while the 
inter-site correlation causes a decrease in width of the majority spin band. There is an 
important difference between these two effects. The situation from Fig. 3 can not lead by 
itself, without the exchange field, to ferromagnetism (see [9,10]), while the majority spin 
bandwidth decrease (see Fig 4) can. However, there is the problem with the magnitude of 
the interaction constants, which can lead to this ordering. There is also a difference in the 
order of approximation necessary to obtain different correlations. The band shape change 
shown in the Fig. 3 can be obtained only in the approximations higher than the I s * order in 
U/D, while the change shown in Fig. 4 is obtained already in the 1 st order approximation 
(Hartree-Fock). 

After inserting Eq. (90) into Eq. (78) one can write for the static ferromagnetic suscep- 
tibility 
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At zero temperature we obtain from Eq. (81) that; It = p°(e F )/b°, where e° F is the 
Fermi level for the system without inter-site interactions and b° is the bandwidth change 
parameter in the paramagnetic state due to the inter-site interactions. This will modify the 
last relations into the following form 



X = 



2/4p°(4) 



(l- Ku -K b )b°-F t F ot p°(e F ) 



(94) 



Comparing the above equation with the equation for susceptibility when there is only 
the Hartree-Fock exchange field 



X 



2/4p°(4) 

1 - F t F ot pO(e° F ) 



(95) 



one can see that the ferromagnetic state will be favored by both positive correlation factors: 
Ku and K b , and by a decrease in the bandwidth b° < 1 due to the inter-site interactions. 

From the zero of the susceptibility denominator (Eq. (94)) and Eqs (80), (23), one can 
calculate the critical on-site exchange interaction. 



per _ 



Kj 



K b )b° 



p{£f) 



zJ -AD 



_d_ 

dm 



sy 



I- a n a -I- a (l-S) 



(96) 



where the kinetic interactions were expressed by the common factor S, see Eq. (5). 
The inter-site correlation factor given by Eq. (92) can be written as 



D 



K h 



-D 

At zero temperature we have P T (s 



db a 

p m =o (e)P T (e) £ g^ de - 



(97) 



K h 



2e° F p°(e F ) db° 



b° 



which produces the following expression 

(98) 



dm 



In further analysis we will use two limiting cases for the on-site Coulomb repulsion - 
weak and strong correlation. 

In the case of weak correlation the Coulomb correlation factor Kjj = 0, and the critical 
on-site exchange interaction depends only on the inter-site interactions through the inter-site 
correlation factor K b . Using I a expressed by Eq. (13) one obtains 



per (l-K b )V 
" P°(4) 



2D 



[l-Sfn 2 



+ (l-5)(l-n) 



(99) 



where the bandwidth change parameter in the paramagnetic state ,6°, is given by 
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1-(1-S)n + (1-S) 



a 2 ( n{2-n) V 2 J + J' - V n{2 - n) 
T~ ^ 4 J t 4 



(100) 



To illustrate the dependence of the critical on-site Stoner field, F cr on electron concen- 
tration (at zero temperature) we will use the above equations and the initial semi-elliptic 
DOS given by 



Mathematically this is the simplest possible DOS which can be used for itinerant ferro- 
magnetism. The rectangular DOS does not have a unique relationship between the Stoner 
field and the magnetization, and the parabolic DOS does not take into account the inter-site 
interactions, since it does not have a finite width and is not centered on the atomic level. 

Electron concentration in the paramagnetic state (at T=0K) is calculated from the simple 
condition 



with expression (101) inserted for the DOS. This will allow to calculate e° F (n) and p°(£%) for 
a given electron occupation n. Finally, we calculate F cr (n) from Eq. (99) with K b estimated 
from Eq. (98), with the help of the derivative |^ calculated from Eq. (22). 

Calculated in this way the dependence of critical on-site Stoner field on electron con- 
centration is shown in the Fig. 5 for the case of weak Coulomb correlation U and various 
parameters of inter-site interactionJ and S, at J' = J and V = 0. Parameters S = 1 and 
J = describe the classic Stoner model. 

We can see from this Figure that the decreasing factor S and increasing interaction, 
J, decrease the on-site Stoner field required to create ferromagnetism, even to zero, at 
some electron concentrations (solid line). Including only the inter-site exchange interaction 
(dashed line) lowers the critical field symmetrically with respect to n = 1. The kinetic 
interactions: t ex and At, described by the hopping inhibiting factor S < 1 (dotted line) 
lower the on-site Stoner field, but now there is no symmetry with respect to the half-filled 
point n — 1. 




(101) 




(102) 
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FIG. 5. The dependence of the critical on-site Stoner field, F cr , on the electron occupation, for 
different values of 5 and J; 5 = 1 and J = (the Stoner model) - dot-dashed line, 5 = 0.6 and 
J = 0.5i - solid line, 5 = 1 and J = 0.5t - dashed line, 5 = 0.6 and J = - dotted line. This is 
the case of the weak Coulomb correlation, U. 



In the case of strong correlation (U » D) the band is split into two sub-bands. For 
the initial semi-elliptic DOS given by Eq. (101) we obtain from Eq. (30) in this limit (see 
[10]) the following densities 

- for the lower Hubbard sub-band 



*D ef \ 



1 - n_ r 



D 



(103) 



for the upper Hubbard sub-band 



/V {£) 



D 



(104) 



with e = e°b°, and D ef = Db° . 

The parameter I a is given by Eqs (16) and (17), which after inserting into Eq. (96) give 
the following values of the critical on-site exchange interaction 

-for the lower Hubbard sub-band 



F „ = (l^Ku z IQl_ zJ _ 2D a-n) 



-for the upper Hubbard sub-band 
_ (1-Ku-K b )b° 



(2-n)' 



zJ -2D(n- 1) 



n 2 + 4 (1 - S) (1 - n) 



(105) 



n- 



(106) 
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where p° (ef) = lim p a (e) with p CT (e) given by Eqs (103), (104) with n_ CT = n/2, and 

m— +0 

£i?(n) calculated from Eq. (102) as Ef{ji) = e F b°. The correlation factor Kjj is calculated 
from Eq. (91), which at T=0K is reduced to the simple condition 



K 



u 



-D, 



dm 



de 



ef 



-D, 



dp-g{e) 
dm 



- F 



de 



dp 



ef 



dm 



-de 



(107) 



in which we will use the DOS given by Eqs (103) and (104). 

We calculate F cr {n) from Eq. (106), with K b given by Eq. (98). The bandwidth b° is 
calculated from Eq. (22) with l a given by the m = limit of Eqs (20) and (21), for the 
lower and upper sub-band, respectively. 

Fig. 6 shows that in the presence of the on-site strong correlation U, the inter-site 
interactions: J, S, decrease the critical Stoner field dramatically. The increase of J causes 
the decrease of F cr for concentrations nearly half-filled and also at small concentrations and 
concentrations close to completely filled. The decrease of S from 1 to causes a drop of 
the critical Stoner field, F cr , especially for small n < 0.5, and for 1 < n < 1.5, where both 
Hubbard sub-bands begin to fill. As we know the strong Coulomb repulsion causes the split 
of the band. This split causes the change in sign of the inter-site correlation factor Kb when 
the Fermi level moves from the lower to the upper Hubbard sub-band. For 5^1 the ratio 
7^ 0. This is why the curves with S < lwill have a discontinuity at half filling, when 



K h \P 

7W) 



the Fermi energy jumps from the lower to the upper sub-band. 




n 



FIG. 6. Dependence of the on-site critical Stoner field, F cr (in the units of half bandwidth), on 
the electron occupation. The curves show the influence of inter-site interactions, J = J', V = 0, 
and of the hopping interactions, At, t ex , represented by the hopping 'inhibiting' factor S on .F cr in 
the presence of strong on-site correlation, U = oo. 



To justify the use in this chapter of the rather academic limit of U » D for F ordering 
we compare on the next graph the three cases; U = 0, U = 3D, and U » D. All the curves 
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were calculated on the base of Eqs (101), (30) and (102). One can see that the difference 
between the curves U = 3D, and U » D is small. This justifies the use in numerical 
calculations (for simplicity) of the strong correlation limit U » D. 




n 

FIG. 7. Dependence of the on-site critical Stoner field, F cr (in the units of half bandwidth), 
on the electron occupation, which shows the influence of on-site correlation, U. The inter-site 
interactions are zero ( J = 0, S = 1) 



Enhancement of magnetic susceptibility 

The susceptibility given by Eq. (93) for T / K or Eq. (94) for T = K is not 
divergent for most of the pure elements, as only few of them are magnetic. But for many 
pure elements the denominator of susceptibility is significantly decreased, which produces 
the experimentally observed increase of susceptibility. This phenomenon was extensively 
studied in the past. From Eq. (94) in the H-F approximation for U, and in the absence of 
the inter-site interactions, when Ky = 0, = 0, and b° = 1 (see e.g. [57]) one can write 



X= 1 F^n^) =XpA ' (1 ° 8) 

c totP \ £ F) 



where the bare Pauli term \p — ^PbP°( £ f) 1S enhanced by the factor A — \j 1 — F[ ot p° (e° F ) 
the Stoner enhancement factor. This enhancement in the case of nonzero correlations is given 
by 



A=l. 



l-Ku-K b )b»-F t F ot p\e» F ) . (109) 



In the past the CPA approximation describing the on-site repulsion U; Ku ^ 0, b° = 1, 
was used to calculate this enhancement (see [9] for the general model of pure elements and 
e.g. [60-62] for susceptibility of disordered binary alloys). Experimental data and theoret- 
ical results from the local density functional method were collected more recently for pure 
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transition elements by [56]. In this chapter the factor A (Eq. (109)) is additionally increased 
by the inter-site interactions; Kj, ^ 0, b° < 1. This effect should be also included in the 
investigation of experimental data on the susceptibility enhancement. 

4.2 NUMERICAL RESULTS FOR MAGNETIZATION AND CURIE TEM- 
PERATURE 

The Curie temperature, at which magnetization vanishes, is calculated from the zero of 
the susceptibility denominator, Eq. (93), attained with temperature 



Quantity F^ ot calculated from this equation at nonzero temperature is larger than F£ t , 
calculated from the same equation (110) at zero temperature, since It has its maximum 
It = p°{s° F )/b at zero temperature. 

In the case of weak Coulomb correlation (U = 0) the Coulomb correlation factor, 
Kjj = 0, hence the Curie temperature depends only on K b and It- 

In the case of strong Coulomb correlation (U » D) according to Eq. (110) we 
have to include the correlation factor K v given by Eq. (91). In this case the band is split 
into two sub-bands. Using the DOS given by Eqs (103) and (104) we calculate numerically 
the correlation factor Ku from Eq. (107) where the self-energy was eliminated 



Alternatively, we can calculate directly the magnetization as; m{T) = n a — n_ a , with 
n± a from Eq. (36) for a given electron occupation; n = n a + ri- a , and then obtain the Curie 
temperature by searching for the Curie point wherem(T c ) — > 0. 

The procedure for calculations is the following. We use Eq. (37) for the electron oc- 
cupation number, which corresponds to a given 3d ferromagnetic element in the Table I. 
We adjust the value of total Stoner field, F[ ot , to create the experimental zero temperature 
magnetization m (T — > K) = n CT (0 K) — ri- a (0 K). Having obtained Ff ot we calculate the 
on-site Stoner field, F, from Eq. (26) at the given assumed inter-site interactions J, S and 
D = 3d half band-width according to [63]. 

Next, we calculate the Curie temperature with the same constants, F, J and S using 
Eq. (38) for the magnetization and searching for the point where m(Tc) — > 0. 

The results in acase of the weak correlation are collected in the Table I. 



\-K v -K h - F?J T = 



(110) 




(111) 



— oo 
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TABLES 



TABLE I. Curie Temperatures and Values of On-site Stoner Field for Ferromagnetic Elements, 
Modified Stoner Model with U = 



Element 


n 


m 


-L>|eVJ 


T C [K] 
F(m)[eV] 








JNo. 5 








t[eV] 


No.l 
5 = 0.6 
J = 0.5t 


No.2 
5 = 0.6 
J = 


No.3 
5 = 1 
J = 0.5t 


No.4 
5 = 1 
J = 




Fe 


1.4 


0.44 


2.8 


2050 


3295 


3980 


4290 


1043 








0.35 


1.18 


3.76 


2.28 


4.85 




Co 


1.65 


0.344 


2.65 


1690 


3300 


3880 


4710 


1388 








0.22 


1.86 


4.46 


3.0 


5.6 




Ni 


1.87 


0.122 


2.35 


620 


870 


1720 


1960 


627 








0.20 


2.33 


4.72 


4.0 


6.4 
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It is very interesting to compare the different results for T c which are shown in this Table. 
Column no. 4, which is the Stoner model for semi-elliptic DOS, shows that all theoretical 
results of the pure Stoner model are much higher than the experimental Curie temperatures 
(column no. 5). This means that the Stoner model, which assumes that the on-site atomic 
field creates ferromagnetism, overestimates, to a large extent, the Curie temperature. The 
necessary on-site Stoner field is also unrealistically high. Perhaps ferromagnetism is the 
result of inter-site forces changing the bandwidth, since column no. 1 is closest to the exper- 
imental results. The larger the component of the on-site field (i.e. S — > 1, J — > 0), the more 
the T c results exceed the experimental values. However, one has to be careful not to overes- 
timate these inter-site interactions since they suppose to be weak. They lower substantially 
the Curie temperature, but the remaining difference between theoretical and experimental 
values is still significant and it has to be attributed to the thermal spin fluctuations, since 
there is an experimental evidence of such fluctuations (see [56,64]) 

These simple calculations would confirm the earlier attempts of understanding the itin- 
erant ferromagnetism as the effect of ordering local moments, whose alignment disappear 
at the Curie temperature. However, the moments themselves exist up to the temperatures 
exceeding the Curie temperature by two or three times (see Mizia [65]). We could think 
about local moments as being created by the Stoner on-site field, but their ordering would 
be driven by the inter-site interactions, which is much weaker than the on-site field. 

In the case of strong correlation the procedure of calculations is the same as in the 
case of the weak correlation. First, we adjust the on-site exchange interaction F(at given 
inter-site parameters J, S) to the magnetic moment at zero temperature. We use the DOS 
given by Eqs (103) and (104). The correlation factor is calculated from Eq. (111). Finally, 
we calculate the Curie temperature from Eq. (110), or from the condition of vanishing 
magnetization, with the same interactions F, J and S. 

Table II shows the results of Curie temperatures obtained for the strong on-site corre- 
lation from the Stoner model (S = 1, J = 0), and for the cases with nonzero kinetic and 
exchange inter-site interactions. It can be seen that the influence of strong on-site correlation 
(U = oo) does not lower the theoretical Curie temperature when compared to the case of no 
on-site correlation (U = 0) (see Table I), although the necessary on-site exchange field, F, 
in the case of strong on-site correlation (U = oo) is lower than without it, U — (compare 
Figs 5 i 6). 
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TABLE II. Curie Temperatures and Values of On-site Stoner Field for Ferromagnetic 3d Ele- 
ments, Modified Stoner Model with U = oo 



Element 


n 


m 


T C [K] 

F(m) [eV] 








T exp^ 








5 = 0.6 


5 = 0.6 


S = l 


S = 1 










J = OM 


J = 


J = 0M 


J = 




Fe 


1.4 


0.44 


2420 
0.41 


4000 
2.45 


3810 
1.03 


4960 
3.06 


1043 


Co 


1.65 


0.344 


1850 
1.35 


3270 
3.65 


3990 
1.82 


5180 
4.12 


1388 


Ni 


1.87 


0.122 


630 
1.94 


910 
4.2 


1770 
2.9 


2020 
5.16 


627 
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Now we want to focus on the temperature dependence of the magnetization m which 
is plotted in Fig. 8. The dependence m(T) was calculated for the strong correlation and 
without the correlation for electron occupation representing iron. In both cases we assumed; 
J = o, S = 1 and J = 0.5t, S = 0.6. 

One can see from this figure that the strong on-site correlation U is not helping fer- 
romagnetism at all, as the curves with U » D lay above the corresponding curves with 
U = 0. 

0.5 1 . . . . . 1 




6000 



T\K] 

FIG. 8. Magnetization dependence on temperature, m(T), without and with strong Coulomb 
correlation, n = 1.4, D = 2.8 eV for different values of parameters J and S; U = 0, J = 0, S = 1- 
solid line, U = oo, J = 0, S = 1- dashed line, U = 0, J = 0.5i, S = 0.6- dotted line, U = oo, 
J = 0.5t, S = 0.6- dot-dashed line. 



5 ANTIFERRO MAGNETISM 



5.1 ONSET OF ANTIFERROMAGNETISM 



We can find the critical value of the total field creating AF from Eq. 
the expression (88) 



(89) using for K x 



K x = J 9 -i^f(e)de= j 9 -^f(e)de 



-oo 
oo 



I /(*) -J Im 



k — oo 



-oo 
9Gg Q (e,k) 
OA 



A^0 



de 



112) 



Differentiating Eq. (64) we can show that 



dG™ (e, k) 



dA 



A^0 



E + Ek E — EkJ ^Ek 



(113) 
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Using this expression in Eq. (112) we obtain the following formula 



f{e k ) - f{-e k ) 
2e k 



1 E- 



i P / \ , P , 
tanh - (e fc - ji) - tanh - (-e fc - fx) 



(114) 



We can calculate now F^ t {n) from the relation; F^ t = 1/K X . The above equation 
although similar to the condition for superconductivity is different in two respects (see [66]). 
First, the summation runs from k minimum in the band to maximum k, giving a positive 
K x and positive F^ t (n). Second, in the denominator we have the dispersion relation e k , 
which is centered on the average energy of the band not around the chemical potential, like 
in the case of superconductivity (see [67]). 

From the expression for the electron numbers n± CT (7 = a, (3) (Eq. (70)) we obtain 
in the mean-field approximation that the antiferromagnetic moment per atom (in Bohr's 
magnetons) is given by the following expression 



m 



= < ~ *%, = ^ E K CT " p r) lf(E k ) - f(-E k )} . 



(115) 



The chemical potential \x is determined from the carrier concentration non the basis of 
the equation, which also comes from Eq. (70) 



n = n 



: + n% = ^J2 (Pt + Pr) [f(E k ) + f(-E k ) 



(116) 



Inserting Eq. (72) into Eq. (115) and using the relation (74) we obtain 

i = -^of^E^[/(^)-/(-^)] • 



(117) 



At the transition from AF to normal state; A — > 0, and the condition (117) takes on the 
form 



1 = 



1 X- /(gfc) ~ /(-£fe) = per 1 V- 1 
tot at 2—/ Or- tot AT 2-~> 



N 



k 



2ei 



N k 4£ k 



, P / \ , P , 
tanh - (e k - fJ>) + tanh - (-e k - n) 



(118) 



where e k = s\b AF , which is the same as the result (114) coming from the static magnetic 
susceptibility. 

For the numerical analysis we will use the A — > limit of Eq. (116) and of Eq. (118) in 
their integral form, 



n 



Deff 

2 J p(e)f(e)de , 

~ D eff 



(119) 
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D, 



D eff 

.,, , , ''I , •^D ~=D^' < 120 > 



—D, 



ff 



FZt I P(e)h tanh § (e - /x) - tanh f (-e - //) 



cfe 



and the semi-elliptic DOS given by Eq. (101) in the case of the weak Coulomb correlation 
U. For the strong Coulomb correlation U we will use split densities given by Eqs (103) 
and (104). The results will show how the critical on-site exchange field depends on the 
electron concentration. Figs 9 and 10 show the dependence of this interaction, F cr , on the 
electron concentration for different values of the inter-site and kinetic interactions, described 
by the parameters; J, S, in the case of weak and strong Coulomb correlation U. Analyzing 
these curves one can see that the inter-site exchange interaction J increases values of on- 
site exchange interaction, F cr , required for AF. The decrease in F cr can be achieved when 
J < 0. This effect is opposite to the case of ferromagnetism, where stabilization of ordering 
was obtained for J > 0. Those two different effects of inter-site exchange interaction on 
magnetism can be roughly understood when one compares the expressions for total field in 
the case of ferromagnetism; F[ ot = F + z[J + t ex n ~ m +2At(l — n)], and antiferromagnetism 
F££ = F — z (J + 2t ex lAF)- We use the word 'roughly' since in calculating the critical 
fields there is included also the factor of the band-width change, which does not appear 
in the effective fields. The result of this simplified approach with respect to the inter-site 
exchange energy, J, is in agreement with the Heisenberg term for the interaction energy of 
localized spins; Si ■ Sj, which points towards ferromagnetism when J > 0, and towards 

antiferromagnetism when J < 0. Analyzing the influence of inter-site interaction J on long 
range ordering one has to remember that the negative value of this interaction stimulates 
the d-wave superconductivity (see e.g. [68]). As a result it will cause the co-existence of SC 
and AF ordering. This is why the t-J model was successful in describing co-existence of the 
singlet SC and antiferromagnetism. Although one has to remember that this model has a 
contradiction in its very origin. This contradiction lies in considering electrons as itinerant 
(the kinetic t-term) and localized (the potential J-term) at the same time. Including in our 
model the kinetic interactions; At and t ex , has decreased the minimum on-site exchange 
interaction necessary for AF. The hopping interaction At decreases F cr by decreasing the 
band- width (see [69]). Positive exchange-hopping interaction t ex also decreases the band- 
width but at the same time it decreases the effective field (see above or Eq. (62)). In effect 
the influence of this interaction on AF is very small. 
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I . . 1 1 1 1 1 1 

0.5 1 1.5 2 

n 

FIG. 9. The dependence of the critical on-site exchange, F cr , on the electron occupation, The 
case of the weak correlation U. The curves for different values of 5 and J are; 5 = 1 and J = 0- 
dot-dashed line, 5 = 0.6 and J = 0.5i - solid line, 5 = 1 and J = 0.5t - dashed line, 5 = 0.6 and 
J = - dotted line. 




n 

FIG. 10. he dependence of the critical on-site exchange, F cr , on the electron occupation. The 
case of the strong correlation U. The curves for different values of S and J are; 5 = 1 and J = 0— 
dot-dashed line, S = 0.6 and J = 0.5t - solid line, 5 = 1 and J = 0.5i - dashed line, 5 = 0.6 and 
J = - dotted line. The double dotted-dashed line is for At = t (ii = or 5 = 0), and t ex = 
(see Eq. (125) below). 



Including into calculations the inter-site interactions; J, J', V, does not shift the minimum 
of the critical field F cr (n), which for the weak correlation is located at n = 1, and for the 
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strong correlation at n — 2/3. The inter-site interactions do not change the character of the 
critical curves, they only lower them by a factor of b AF from Eq. (39) or (40). Remembering 
that £ = e°b AF this can be proved analytically 



d. 



eff 



K x (n) = - J p(e) /(£) -f ( - £) & 

-Deff _ ( 12 1) 

= ■ j p o (e0) SM^t. = ■ Kl („) 



b - D 26 

Since F% t = 1/K X , and F cr = l/K° x we obtain that 

=F cr -z(J + 2t ex I AF ) = b AF Fq T , (122) 

where F cr and Fq t are the critical on-site exchange fields for AF with and without the inter- 
site interactions, respectively. 

Including the kinetic interactions, S < 1, leaves the minima of critical curves where they 
are as it only changes the factor b AF in Eq. (122). In this chapter we assumed dependence 
of the hopping energy on the occupation expressed by the equation 

t\j = t - At(ni^ a + fij-a) + 2t ex n i - a n j - (J . (123) 

This relation together with the other inter-site interactions; J, J', V, expressed in the 
generalized Hartree-Fock approximation brought the Eq. (39) for b AF with Iaf given by Eq. 
(20) or (21) in the case of U » D. The only exception from the scaling rule for critical 
on-site exchange interaction in the case of U » D is the situation when b AF (n = 1) = 0. 
This would cause the critical curve to drop at n = 1, what would make antiferromagnetism 
possible at this concentration. In the itinerant band model this is the situation when there 
is localization at half-filling. As a result we obtain at this concentration possible AF, and 
since the band will expand very rapidly with the occupation n departing from one there 
will be also a strong driving force towards SC (see [70]). In such circumstances we will 
have the co-existence or rather competition at n — 1 between AF and SC. Unfortunately, 
sincei ' af (n — 1) = 0, the relation (39) givesb AF (n — 1) ^ 0. In our previous paper (see 
[69]) we have the zero bandwidth at n — 1, but we have neglected the t ex term in Eq. (3) or 
(123) above. This is equivalent to the use of Hirsch's simplified linear approach (see [47,71]) 

= t - At(n^ a + n^ a ) , (124) 

instead of Eq. (123). From this relation it follows that 

bAF ^ 1 _M n _2_J ± J>-V lAF (125) 

At At = t or ti = this equation givesfe^ (n = 1) = 0, and we have the possibility of 
AF ordering, see the double dotted-dashed line in Fig. 10. AF can occur at n — 1 only 
when t\ = what means that there is no hopping in the presence of other electrons with 
opposite spin, and additionally when t ex = 0, which leads to the strange condition that 
t 2 = —t. Alternatively we can obtain the Eq. (124) (which leads to the AF at n — 1) from 
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the basic Eq. (2) in a first approximation by neglecting all the operator products of the 
type; h^n^ in Eq. (2). 

5.2 NUMERICAL RESULTS FOR MAGNETIZATION AND NEEL'S TEM- 
PERATURE 

The analysis of Eqs (116) and (117) in their integral form will give the magnetization 
dependence on the temperature for different values of the on-site and inter-site interactions. 

D 

n= J p(e) [f(E) + f(-E)]de , (126) 

—D 



D 



Faf / p{e) 



-D 



[f(E) - f(—E)} 
2E 



de 



(127) 



where 



E = Ve 2 + A 2 



A 



F AF- 

r tot 2 



711 



7 AF 
tot 



F-z{J + 2t ex I AF ) 



(128) 



From the above equations we find numerically the sub-lattice magnetization, and the 
Neel's temperature, which is the temperature where this magnetization drops to zero. Fig. 
11 presents the dependence of Neel's temperature, Tjy, on electron concentration for different 
values of the inter-site and kinetic interactions for the weak Coulomb correlation. All the 
curves are for the same on-site exchange field. The curves with nonzero inter-site exchange 
interaction J have a lower Neel's temperature, because this interaction decreases the effective 
exchange interaction; F^f, F^f = F — z (J + 2t ex lAF)- This effect is stronger than the 
decrease of the bandwidth due to J (see Eq. (39)), which increases T N . The assisted 
hopping interaction At = t(l — S), and kinetic-exchange interaction t ex = |(1 
the Neel temperature by decreasing the bandwidth (Eq. (39)). 
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FIG. 11. Dependence of the Neel's temperature (T/v) on electron concentration for D = 0.5 eV, 
F = 0.34 eV and for different values of 5 and Jin the case of the weak Coulomb correlation; 5 = 1 
and J = - dot-dashed line, 5 = 0.7 and J = 0.2i - solid line, 5 = 1 and J = 0.2t - dashed line, 
5 = 0.7 and J = - dotted line. 



Figure 12 shows the magnetization versus temperature for electron occupation n = 0.95 
and for various values of the inter-site and kinetic interactions. Again, all the curves are 
calculated for the same internal exchange field F = 0.34 eV. Similarly to the calculations 
of Neel temperature (in Fig. 11) at constant value of the on-site exchange field, the highest 
magnetization is obtained for parameter S < 1, which means that magnetization is enhanced 
when electron hopping is inhibited in the presence of other electrons. 




T/T N 

FIG. 12. Magnetization versus temperature for band filling n = 0.96, D = 0.5 eV, F = 0.34 eV 
and for different values of 5 and J; 5 = 1 and J = - dot-dashed line, 5 = 0.7 and J = 0.2t - 
solid line, 5 = 1 and J = 0.2i - dashed line, 5 = 0.7 and J = - dotted line. This is the case of 
the weak Coulomb correlation. 



TABLE III. Neel's Temperatures and Values of On-site Stoner Field for Antiferro- 
magnetic Elements, U = 



Element n 



D [eVl Tn[K] 

U [eWl F AF (m) [eV] 



No. 5 
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No.l 


No.2 


No.3 


No. 4 








S = 0.6 


S = 0.6 


S= 1 


S = 1 








J — u.ot 


t — n 
J — u 


j — yj.oz 


t — n 




Cr 1.08 


0.08 3.5 


520 


1430 


1530 


2440 


311 [72] 






2.37 


1.47 


3.17 


2.27 




Mn 1.24 


0.48 2.8 


2110 


5840 


7090 


10820 


540 [56] 






2.41 


2.65 


4.47 


4.71 
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Table III shows the Neel's temperature calculated for the 3d antiferromagnetic elements 
(Cr and Mn) in the case of the weak correlation. First the on-site exchange interaction, F, at 
a given assumed inter-site interactions J and S, was found from Eq. (43) after fitting F^f to 
the experimental magnetic moment at T = 0. Next, using the same F the Neel's temperature 
was calculated. The results show that the inter-site, J, and kinetic, S, interactions decreased 
the Neel's temperature (column no. 1). The difference between experimental and theoretical 
Neel's temperature is even larger than in the case of ferromagnetism (confer Table I). The 
large difference has to be attributed, as in the case of F, to the spin waves or spin fluctuations 
(see [56,64,72]). 

6. CONCLUSION 

We will summarize now the results for ferromagnetism and antiferromagnetism obtained 
in this article. Our analysis of the Hubbard model was focused on the influence of the 
on-site and inter-site Coulomb correlation and the kinetic interactions on magnetism. We 
analyzed magnetic susceptibility, the values of critical exchange interaction, the transition 
temperature and the magnetization. 

In analyzing the ferromagnetic ordering we have found three driving forces contributing 
to this ordering. 

(i) spin dependent band shift (Stoner shift) coming from the on-site F and inter-site 
exchange interactions J, V and the kinetic interactions At, t ex (represented by the parameter 
S), 

(ii) spin dependent change of the band width and band capacity depending on the electron 
concentrations due to the on-site Coulomb repulsion described by the on-site correlation 
factor Kjj. Increased capacity of the majority spin sub-band drives the magnetism. 

(iii) spin dependent change of the band width and band shape depending on the electron 
concentrations due to the inter-site interactions described by the inter-site correlation factor 
Kb . Narrowing of the majority spin sub-band drives the magnetism. 

The first effect is the classic effect known since the approach of Weiss [73] and Slater 
[74] to magnetism. This is lowering of the potential exchange energy during the transition 
from the paramagnetic to ferromagnetic state. During this transition there is an increase in 
the kinetic energy. The balance of the sum of these two energies decides whether the actual 
transition takes place. This balance leads to the existence of the critical values for different 
interactions, above which the transition takes place. In its original form it gives us the Stoner 
condition for ferromagnetism; F^,p(eF) > 1- In the new development, described in points 
two and three above, the increase in the kinetic energy is moderated by the on-site and 
inter-site correlation effects. These correlations will also change (decrease) different critical 
interactions coming from the total energy balance. The inter-site and on-site correlations 
are coming into the critical condition throughout the correlation factors, which appear in the 
denominator of the susceptibility (Eq. (84)) and as the coefficient of second order term in 
the free energy expansion over magnetization (Eq. (85)), which is proportional to the inverse 
of susceptibility. The inter-site and on-site correlations reduce significantly the denominator 
of the static magnetic susceptibility (Eq. (84)). For this reason the magnetic susceptibility 
is enhanced for the materials with strong electron correlation (Eq. (109)). The zeros of the 
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denominator give the critical values of the total interaction, which is greatly reduced with 
respect to the classic Stoner model. 

Analysis of the inter-site correlation has shown the great decrease of the Curie tempera- 
ture towards the experimental values (see Table I i II). First, the direct calculations within 
the original Stoner model (see Section 4.2) have pointed out that the Curie temperature, 
after fitting the Stoner shift to the experimental magnetic moment at zero temperature, is 
much too high (see Table I). Next, using the modified H-F approximation for the inter-site 
interactions, which is equivalent to introducing the inter-site correlation in the first order 
approximation, we arrived at much lower Curie temperatures. Considering the simplicity 
of the model, and the additional presence of spin waves which are out of the scope of this 
chapter, the results were close enough to the experimental data (see Table I). Apparently 
the inter-site interactions are 'softer' and decrease faster with the temperature than the 
on-site Stoner field used in the original Stoner model (see also [75]). Adding up the on-site 
strong correlation U to the Stoner model do not improve the values of the theoretical Curie 
temperature (see Table II). This fact can be understood better after examining Fig. 6 and 
Fig. 7, where the critical field (initializing magnetization) has dropped down at half-filling 
but not at the end of the band, where the 3d elements are located. 

As already established in the past the electron on-site correlation can help in creating 
antiferromagnetism (AF) at half-filled band, where the antiferromagnetic 3d elements 
(Cr, Mn) are located, by dropping the critical field for AF to zero (see Fig. 9) in the case of 
the weak on-site correlation. The strong U, which gives the correct results for the cohesion 
energy in the middle of the 3d row (see [63,76,77]), shifts the zeros of the critical field for 
AF to the maxima of the two split sub-bands, which are located at n = 2/3, and n = 4/3 
(see Fig. 10). Unfortunately, this does not agree with the experimental evidence for AF, 
which is present around n ~ 1. The introduction of inter-site interactions in the modified 
H-F approximation does not shift these minima. It only lowers the critical curves (see Fig. 
10). The inter-site interactions (J,J',V) do not reduce the bandwidth at n — 1 (which 
would lower the critical curves at this concentration), in the case of the strong Coulomb 
correlation, since Iaf (n — 1) — (see Eqs (20) and (21)). 

There is one exception to this rule. We can obtain AF at half-filling at large U if the 
bandwidth at n — 1 goes to zero. This can be achieved in the linear approximation, which 
was described in §5.1. The bandwidth was reduced to zero by assuming that the hopping 
integral in the presence of another electron is forbidden at concentration n — 1; t\ — 0, and 
additionally that t ex = (see Eq. (125)). 

In conclusion the antiferromagnetism at half-filling in the high correlation case can appear 
when the bandwidth goes to zero at this concentration. It can go to zero either as: 1 — n^, 
due to the presence of the assisted hopping interaction or due to some other interaction, 
whose strength is proportional to n. 

This situation can describe the high temperature superconducting cuprate YBaCuO, 
where at the half-filling and the strong Coulomb correlation there is, initially, the anti- 
ferromagnetic order. The superconductivity will appear upon doping only away from the 
half-filling at the electron concentration n ~ 0.95 — 0.8 (see e.g. [78]). 

The temperature dependence of AF was analyzed numerically in Section 5.2. Again, as 
in the case of ferromagnetism the direct calculations within the original Stoner model have 
pointed out that the Neel's temperature, after fitting the Stoner shift to the experimental 
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magnetic moment at zero temperature, is much too high (see Table III). Next, using the 
modified H-F approximation for the inter-site interaction, we arrived at much lower Curie 
temperatures. They are not as close to the experimental data as in the case of ferromag- 
netism. The reason for this could be the existence of non linear modes of excitations in 
Chromium and Magnesium (see e.g. [72]). Nevertheless, as already mentioned above, the 
model is very simple, the details of the realistic DOS as well as the magnetization decrease 
through the spin waves excitation or moment fluctuations could be included, which would 
bring the theoretical results to complete agreement with the experimental data. 



40 



REFERENCES 



[1] Hubbard J. Proc. R. Soc. A 1963, 276, 238. 
[2] Gutzwiller M.C. Phys. Rev. Lett. 1963, 10, 159. 
[3] Kanamori J. Prog. Theor. Phys. 1963, 30, 275. 
[4] Stoner E.C. Proc. Roy. Soc. A 1936, 154, 656. 
[5] Stoner E.C. Proc. Roy. Soc. A 1938, 165, 372. 
[6] Wohlfarth E.P. Philos. Mag. 1951, 42, 374. 
[7] Hubbard J. Proc. Roy. Soc. 1964, A281, 401. 

[8] Velicky B.; Kirkpatrick S.; Ehrenreich H. Phys. Rev. 1968, 175, 747. 

[9] Fukuyama H.; Ehrenreich H. Phys. Rev. B 1973, 7, 3266. 
[10] Mizia J. phys. stat. sol. (b) 1976, 74, 461. 
[11] Nolting W.; Borgie W. Phys. Rev. B 1989, 39, 6962. 
[12] Herrmann T.; Nolting W. J. Magn. Magn. Mater. 1997, 170, 253. 
[13] Herrmann T.; Nolting W. Phys. Rev. B 1996, 53, 10579. 
[14] Hirsch J.E. Phys. Rev. B 1985, 31, 4403. 

[15] Pruschke T.; Jarrell M.; Freericks J.K. Adv. Phys. 1995, 44, 187. 

[16] Georges A.; Kotliar G.; Krauth W.; Rozenberg M.J. Rev. Mod. Phys. 1996, 68, 13. 

[17] Metzner W.; Vollhardt D. Phys. Rev. Lett. 1989, 62, 324. 

[18] Blankenbecler R.; Scalapino D.J.; Sugar R.L. Phys. Rev. D 1981, 24, 2278. 

[19] Hirsch J.E.; Fay R.M. Phys. Rev. Lett. 1986, 56, 2521. 

[20] Jarrell M. Phys. Rev. Lett. 1992, 69, 168. 

[21] Wahle J.; Bliimer N.; Schlipf J.; Held K.; Vollhardt D. Phys. Rev. B 1998, 58, 12749. 

[22] Ulmke M. Eur. Phys. J. B 1998, 1, 301. 

[23] Meyer D.; Nolting W. Phys. Rev. B 2000, 62, 5657. 

[24] Wegner T.; Potthoff M.; Nolting W. Phys. Rev. B 1998, 57, 6211. 

[25] Byczuk K.; Ulmke M.; Vollhardt D. Phys. Rev. Lett. 2003, 90, 196403. 

[26] Amadon J.C.; Hirsch J.E. Phys. Rev. B 1996, 54, 6364. 

[27] Hirsch J.E. Phys. Rev. B 1999, 59, 6256. 

[28] Hirsch J.E. Phys. Rev. B 1989, 40, 2354. 

[29] Hirsch J.E. Phys. Rev. B 1989, 40, 9061. 

[30] Hirsch J.E. Phys. Rev. B 1991, 43, 705. 

[31] Kollar M.; Strack R.; Vollhardt D. Phys. Rev. B 1996, 53, 9225. 

[32] Strack R.; Vollhardt D. Phys. Rev. Lett. 1994, 72, 3425. 

[33] Gorski G.; Mizia J. Physica B 2004, 344, 231. 

[34] Gorski G.; Mizia J.; Kucab K. Physica B 2003, 325, 106. 

[35] Vollhardt D.; Bliimer N.; Held K.; Kollar M.; Schlipf J.; Ulmke M.;Wahle J. Adv. Solid 

State Phys. 1999, 38, 383. 
[36] Kollar M.; Vollhardt D.Phys. Rev. B 2001, 63, 045107. 

[37] Barreteau C; Desjonqueres M.C; Ole A.M.; Spanjaard D. Phys. Rev. B 2004,69, 
064432. 

[38] , Saxena S.S.; Agarwal P.; Ahilan K.; Grosche F.M.; Haselwimmer R.K.W.; Steiner 
M.J.; Pugh E.; Walker I.R.; Julian S.R.; Monthoux P.; Lonzarich G.G.; Huxley A.; 
Sheikin L; Braithwaite D.; Flouquet J. Nature 2000, 406, 587. 

[39] Aoki D.; Huxley A.; Ressouche E.; Braithwaite D.; Flouquet J.; Brison J-P.; Lhotel E.; 
Paulsen C. Nature 2001, 413, 613. 



41 



[40] Pfleiderer C; Uhlarz M.; Hayden S.M.; Vollmer R.; Lohneysen H.V.; Bernhoeft N.R.; 

Lonzarich G.G. Nature 2001, 412, 58. 
[41] Spalek J.; Wrobel R; Wojcik W. Physica C2003, 387, 1. 
[42] Mizia J.; Gorski G.; Kucab K. phys. stat. sol. (b) 2002,229, 1221. 
[43] Kyung B. Phys. Rev. B 2000, 62, 9083. 

[44] Micnas R.; Ranninger J.; Robaszkiewicz S. Rev. Mod. Phys. 1990, 62, 113. 

[45] Campbell D.K.; Gammel J.T.; Loh Jr E.Y. Phys. Rev. B 1988, 38, 12043. 

[46] Campbell D.K.; Gammel J.T.; Loh Jr E.Y. Phys. Rev. B 1990, 42, 475. 

[47] Hirsch J.E. Phys. Rev. B 1993, 59, 3327. 

[48] Hirsch J.E. Phys. Rev. B 2002, 65, 184502. 

[49] Gunnarsson O.; Christensen N.E. Phys. Rev. B 1990, 42, 2363. 

[50] Foglio M.; Falicov L. Phys. Rev. B 1979, 20, 4554. 

[51] Aligia A. A.; Gagliano E.; Arrachea L.; Hallberg K. Eur. Phys. J. B 1998, 5, 371. 

[52] Soven P. Phys. Rev. 1967, 156, 809. 

[53] Plischke M.; Mattis D.Phys. Rev. B 1973, 7, 2430. 

[54] Brouers F.; Giner J.; Van der Rest J. J. Phys. F 1974, 4, 214. 

[55] Mizia J. phys. stat. sol. (b) 1977, 84, 449. 

[56] Kiibler J. Theory of Itinerant Electron Magnetism; International Series of Monographs 

in Physics; Clarendon Press, Oxford, GB, 2000. 
[57] Mohn P. Magnetism in the Solid State; Spinger-Verlag, Berlin Heidelberg, , 2003. 
[58] Kittel Ch. Introduction to Solid State Physics; John Wiley and Sons, New York, Chices- 

ter, Brisbane, Toronto, Singapure 1996. 
[59] Morrish A.H. The Physical Principles of Magnetism; John Wiley and Sons Inc., New 

York, US, 1965. 

[60] Levin K.; Bass R.; Bennemann K.H. Phys. Rev. B 1972, 6, 1865. 
[61] Fukuyama H. Phys. Rev. B 1972, 5, 2872. 
[62] Hasegawa H.; Kanamori J. J. Phys. Soc. Jap. 1971, 31, 382. 
[63] Kajzar F.; Mizia J. J. Phys. F 1977, 7, 1115. 

[64] Moriya T. Spin Fluctuations in Itinerant Electron Magnetism; Springer Verlag, Berlin, 
DE, 1985. 

[65] Mizia J. J. Phys. F: Met. Phys. 1982, 12, 3053. 

[66] Mizia J.; Romanowski S.J. phys. stat. sol. (b) 1994, 186, 225. 

[67] Bardeen J.; Cooper L.N.; Schrieffer J.R. Phys. Rev. 1957, 108, 1175. 

[68] Micnas R.; Ranninger J.; Robaszkiewicz S. Phys. Rev. B 1989, 39, 11653. 

[69] Gorski G.; Mizia J.; Kucab K. Physica B 2003, 336, 308. 

[70] Hirsch J.E.; Marsiglio F. Phys. Rev. B 1989, 39, 11515. 

[71] Hirsch J.E. Phys. Rev. B 2000, 62, 14487. 

[72] Fawcett E.Rev. Mod. Phys. 1988, 60, 209. 

[73] Weiss P. J. de Physique 1907, 6, 661. 

[74] Slater J.C. Phys. Rev. 1936, 49, 537. 

[75] Hubbard J. Phys. Rev. B 1979, 19, 2626. 

[76] Friedel J.J. Phys. 71978, 39, 651. 

[77] Kajzar F.; Friedel J.J. Phys. 1 1978, 39, 397. 

[78] Tallon J.L.; Bernhard C; Shaked H.; Hitterman R.L.; Jorgensen J.D. Phys. Rev. B 
1995, 51, 12911. 



42 



